Oscillation bands of condensates on a ring: Beyond the mean field theory 
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Abstract: The Hamiltonian of a iV-boson system confined on a ring with zero spin and repulsive 
interaction is diagonalized. The excitation of a pair of p- wave-particles rotating reversely appears to 
be a basic mode. The fluctuation of many of these excited pairs provides a mechanism of oscillation, 
the states can be thereby classified into oscillation bands. The particle correlation is studied 
intuitively via the two-body densities. Bose-clustering originating from the symmetrization of wave 
functions is found, which leads to the appearance of 1-, 2-, and 3-cluster structures. The motion is 
divided into being collective and relative, this leads to the establishment of a relation between the 
very high vortex states and the low-lying states. 



After the experimental realization of the Bose-Einstein 
condensation 1 , various condensates confined under dif- 
ferent circumstances have been extensively studied the- 
oretically and experimentally. Mostly, the condensates 
are considered to be confined in a harmonic trap. Con- 
densates trapped by periodic potential have also been 
studied due to the appearance of optical lattices. 2 It 
is believed that the appearance of condensates confined 
in particular geometries is possible. Experimentally, the 
particle interactions can now be tuned from very weak 
to very strong, 3-8 it implies that the particle correlation 
may become important. Theoretically, to respond, go- 
ing beyond the mean field Gross-Pitaevskii (GP) theory 
is desirable, and the condensates confined in particular 
geometries are also deserved to be considered. 

Along this line, in addition to the ground state, the 
yrast states have been studied both analytically and 
numerically. 9-17 The condensation on a ring has also 
been studied recently. 12 The present paper is also ded- 
icated to the N— boson systems confined on a ring with 
weak interaction, its scope is broader and covers the 
whole low-lying spectra. A similar system has been in- 
vestigated analytically by Lieb and Liniger 16 ' 17 . How- 
ever, the emphasis of their papers is different from the 
present one, which is placed on analyzing the structures 
of the excited states to find out their distinctions and sim- 
ilarities, and to find out the modes of excitation. Based 
on the analysis, an effort is made to classify the ex- 
cited states. Traditionally, the particle correlation and 
its effect on the geometry of N— boson systems is a topic 
scarcely studied if N is large. In this paper, the corre- 
lation is studied intuitively so as the geometric features 
inherent in the excited states can be understood. Tradi- 
tionally, a separation between the collective and internal 
motions is seldom to be considered if N is large. In this 
paper such a separation is made and leads to the estab- 
lishment of a relation between the vortex states and the 
low-lying states. 



It is assumed that the N identical bosons confined on 
a ring have mass to, spin zero, and square-barrier inter- 
action. The ring has a radius R, N is given at 100, 20 
and 10000. Let G = h 2 /{2mR 2 ) be the unit of energy. 
The Hamiltonian then reads 
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(1) 



where 0i is the azimuthal angle of the i-th boson. r%J 



Vii 

V if \6j — 6i\ < 9 r an g e , or = otherwise. Let 4>k — 
e lke /VStt be a single particle state, — fc max < k < k max 
is assumed. The iV— body normalized basis functions 
in Fock-representation are \a) = |?i-fc max ,' ■ • "fc max ), 
where rij is the number of bosons in <j)j , y\ rij = N, 
and rijj — L, the total angular momentum. Then, 
H is diagonalized in the space spanned by \a), the low- 
lying spectrum together with the eigen- wave- functions, 
each is a linear combination of \a), are thereby obtained. 
Let K a = J2j n jj 2 be the total kinetic energy of an \a) 
state. Evidently, those \a) with a large K a are negligi- 
ble for low-lying states. Therefore, one more constraint 
K a < ivmax is further added to control the number of 
| a). In this procedure, the crucial point is the calcula- 
tion of the matrix elements of H. This can be realized 
by using the fractional parentage coefficients 18 (refer to 
eq.(6) below). Numerical results are reported as follows. 

This paper concerns only the cases with weak interac- 
tion. Firstly, let V a = l,9 range = 0.025, and N = 100. 
This is corresponding to 7 = 0.00157 , where 7 is in- 
troduced by Lieb and Liniger to measure the strength 
of interaction, 16,17 this is shown later. When fc max and 
if max are given at a number of values, the associated 
eigen-energies Ej of the first, fifteenth, and sixteenth L = 
eigen-states are listed in Table I. When (A: max , K max ) 
is changed from (3, 50) to (5,60), the total number of \a) 
is changed from 2167 to 8890. Table I demonstrates that 
the great increase of basis functions does not lead to a 
remarkable decrease of eigen-energies. Thus the conver- 
gency is qualitatively satisfying even for the higher states. 
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TABLE I: Eigen-energies E 3 (the unit is G) of the L — 
states. The first row is (fc ma x, /fmax), the first column is the 
serial number of states j. V = I, Orange = 0.025, and N = 
100 are given. 





(3,50) 


(4,50) 


(4,60) 


(5,60) 


1 


39.109 


39.090 


39.090 


39.078 


15 


53.645 


53.616 


53.616 


53.613 


16 


54.822 


54.800 


54.800 


54.790 



In the following the choice /c max = 4 and K max = 50 
are adopted, this limitation leads to a 3254-dimensional 
space. Thereby the resultant data have at least three 
effective figures, this is sufficient for our qualitative pur- 
pose. 



N=100, L=o 



FIG. 

G = 



1: The spectrum of L = states, the unit of energy is 



H 2 /(2mR 2 ). N = 100, V = 1, and 6 r 



0.025 are 



assumed, they are the same for Fig.l to Fig. 5. The levels in a 
column constitute an oscillation band, the levels in bold line 
are doubly degenerate. 

The low- lying spectrum is given in Fig.l, where the 
lowest fourteen levels are included. Twelve of them can 
be ascribed into three bands, in each band the levels 
are distributed cquidistantly, this is a strong signal of 



(L) 
Z,i 



harmonic-like oscillations. From now on the labels 'F 

and E( L ' Z ' 1 ^ are used to denote the wave function and en- 
ergy of the i-th state of the Z-th band (Z=I, II, III,- • •). 

It turns out that the excitation of a pair of particles 
both in p-wave but rotating reversely, namely, one par- 
ticle in <j>i while the other one in (p-i, is a basic mode, 
the pair is called a basic pair in the follows. A number of 
such basic pairs might be excited. When 2j particles are 
in basic pairs while the remaining N — 2j particles are 
in (j>Q, the associated \a) is written as \P^). For all the 
states of the /—band, we found \e4°] is mainly a linear 
combination of \P^) together with a small component 
denoted by Aj^, i.e., 



(2) 



where Aj ; , is very small as shown in Table II, while the 
coefficients cj L ' z '*) arise from the diagonalization. Thus 



the basic structure of the /—band is just a fluctuation of 
many of the basic pairs. 

TABLE II: The weights of A z ,i of the bands with L = 



/—band 



//—band 



///-band 



0.009 
0.012 
0.021 
0.035 
0.055 
0.079 



0.017 
0.030 
0.061 
0.028 
0.035 
0.106 



0.040 
0.056 
0.088 



For lower states, Cj°' Z ' % ^ would be very small if j 
is larger, e.g., for the ground state, Cq ' 1 ' 1 " 1 = 0.968 
and Cj°>2^ ~ 0j if implies that the excitation of many 
pairs is not probable. It also implies that the ground state 
wave function obtained via mean-field theory might be a 
good approximation. However, for higher states, many 
pairs would be excited. E.g., for the third state of the 
/-band, Cf' 7,3) = 0.052, 0.406, 0.681, -0.528, and 0.245 
when j is from to 4, it implies a stronger fluctuation. 

When a \a) has not only 2j particles in the basic pairs, 
but also m particles in (fik, while the remaining particles 
in (f>o, then it is denoted as \{k) m P^) (where k = ±1 
are allowed) Similarly, we can define |(fci) mi (fc 2 ) m2 P^), 
and so on. For all the states of the //—band, we found 



vr/o) __ r (0,II,i) 1 r 



±|(-2) 1 (1) 2 P«) 



(3) 



in. i 



where both the + and — signs lead to the same energy, 
thus the level is two-fold degenerate. Again, all the A// ; j 
are very small as shown in Table II, thus the fluctua- 
tion of basic pairs is again the basic structure. However, 
the //—band is characterized by having the additional 
3-particle-excitation (one in d-wave and two in p-wave). 
For all the states of the third band, we found 



(0) 
III, i 



{0,1 1 1, i) 



|(2) 1 (-2) 1 P«)+A /7 ^ (4) 



Thus, the ///—band contains, in addition to the fluc- 
tuation of basic pairs, a more energetic pair with each 
particle in d-wave. It was found that the spacing 
£j(o,z,i+i) _ £{o,z,i) ingidg a n ^ ne bands are nearly the 

same, they are ^3.15. This arises because they have the 
same mechanism of oscillation, namely, the fluctuation of 
basic pairs. 

When the energy goes higher, more oscillation bands 
can be found. The two extra levels in Fig.l at the right 
are the band-heads of higher bands. 

Incidentally, the band-heads of the above three 
bands are dominated by |(2) 1 (-f) 2 p(°)) ± 

|(-2) 1 (1) 2 P(°)) and |(2) 1 '(-2) 1 P(°)), respectively, and 
their kinetic energies K a = 0,6, and 8. Among all the 
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basis functions with L = and without basic pairs, these 
three are the lowest three. This explains why the band- 
heads are dominated by them. Once a band-head is fixed, 
the corresponding oscillation band would grow up via the 
fluctuation of basic pairs. 

The particle correlations can be seen intuitively by ob- 
serving the two-body densities 

P2{6i,0 2 ) = jd6z---de N *<g**<g (5) 

Similar to the calculation of the matrix elements of 
interaction, the above integration can be performed in 
coordinate space by extracting the particles 1 and 2 
from | a) by using the fractional parentage coefficients 18 , 
namely, 

\a) = E fe Vn k (n k - 1)/N(N - i)<fe(l)<fc(2)jq fc ) 
+E ka , kb V n k a n kb /N(N - l)0 ka (l)<t> kb (2)\a kakb ) 

(k a ^k b ) 

(6) 

where \a k ) is different from \a) by replacing n k with n k — 
2, \ a k a k b ) is different from \a) by replacing n ka and n kb 
with n ka — 1 and n kb — 1, respectively. 
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FIG. 2: p 2 as functions of 2 for the 1(a), 11(b), and III (c) 
bands of L — states, #i = is given. The labels i of the 
states j are marked by the curves. 

P2 gives the spatial correlation between any pair of 
particles as shown in Fig. 2. For the ground state ^^j, 
P2 is flat implying that the correlation is weak. However, 
it is a little larger when the two particles are opposite to 
each other (6\ = and 62 — n). It implies the existence 
of a weak correlation which is entirely ignored by the 
mean field theory. Thus, even the interaction adopted is 
weak and even for the ground state, there is still a small 
revision to the mean field theory. For higher states of the 
/—band, the fluctuation of basic pairs becomes stronger. 
Due to the fluctuation, the particles tend to be close to 
each other to form a single cluster. This tendency is 
clearly shown in Fig. 2a. 

For the first state of the III — band, ^jjj 1 has two 
peaks in pi implying a 2-cluster structure. It arises from 



the two d-wave paticles inherent in the band. The fea- 
ture of ^f} x is lying between and \t4jj v For all 
higher states of every band, due to the strong fluctuation 
of basic pairs, all the particles tend to be close to each 
other as shown in 2b and 2c. 

To understand the physics why the particles tend to be 
close to each other, let us study the most important basis 
state \P(ri). By inserting \P^) into eq.(5) to replace 
^ { zl and by using (6), P2 reads 

P2(6i,e 2 ) = w^N-v WN-V-jiW-Gj) 
+4j(N - 2j)(l + cos(6 1 - 6 2 )) + 4j 2 cos 2 (0! - 6 2 )\ 

(7) 

Where there are four terms at the right, the non- 
uniformity arises from the third and fourth terms. The 
third term causes the particles to be close to each other 
to form a single cluster, while the fourth term causes 
the two-cluster clustering. When j is small, the fourth 
term can be neglected, and the particles tend to form 
a single cluster. However, when j w N/2 , the third 
term can be neglected, and the particles tend to form 
two clusters. It is noted that, if the symmetrization were 
dropped, the density contributed by \P^) would be uni- 
form. The appearance of the clustering originates from 
the symmetrization of the bosonic wave functions, there- 
fore it can be called as bose-clustering. 

For L = 1 states, the lowest energy E^ 1 ' 1 ^ is higher 
thanE 1 ^- 7 ' 1 ) by 1.606, but lower than E^' 1 ^ . Thus $ 
is the true first excited state of the system. A number 
of oscillation bands exist as well, the wave functions of 
the lowest six bands are found as 

*$ = EiCf ,J,<) 1(1)^) + Aj^ 

*//< - Ej Cf- H ' l) |(2) 1 (-1) 1 P (J) > + Ajj.i 

= £j Cf ,m><) |(-2) 1 (1) 3 P^)) + Am, 
= Ej Cf'' V,<) |(2) 1 (-2) 1 (1) 1 J P°' ) > + A/v.i 
^ { v\ = T, j Cf' V ' i) |(3) 1 (-l) 2 i > ^> + A v , i 

^vl i = Ej cf yi ' l) \(3)H-2ypU)) + A VIti 

(8) 

Where the weights of all the A z , t < 0.1 if i < 4. Thus, 
just as the above L = case, all the bands have the 
common fluctuation of basic pairs, but each band has 
a specific additional few-particle excitation. The ener- 
gies of the band-heads from I to VI are 40.70, 45.42, 
48.58, 50.14, 52.04, and 53.58 respectively. Further- 
more, the spacing ~3.15 found above is found again for 
all these bands due to having the same mechanism of os- 
cillation. The /—band is similar to the above /—band 
with L = but having an additional single p-wave ex- 
citation, the P2 of them arc one-one similar. Similarly, 
the p2 of the IV— band is one-one similar to those of 
the above ///— band with L = 0. The p2 of the // 
and ///—bands are both similar to those of the above 
//—band with L = 0. However, the V and VI— bands 
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are special due to containing the f-wave excitation, the pi 
of their band-heads exhibit a 3-cluster structure as shown 
in Fig. 3. When the energy goes even higher, more higher 
oscillation bands will appear. For the above six L = 1 
bands, their band-heads are dominated by the |a) with 
K a = 1, 5, 7, 9, 11, and 13. Obviously, a higher K a leads 
to a higher band. 




120 18 



240 300 



FIG. 3: p2 for selected L — 1 states. 8\ = 0, the (Z, i) labels 
are marked by the curves. 

In general, all the low-lying states can be classified into 
oscillation bands. For all the lower bands disregarding 
L, it was found that each band-head is dominated by 
a basis function containing a specific few-particle exci- 
tation but not containing any basic pairs. The energy 
order of the bands is determined by the magnitudes of 
K a associated with the dominant basis function \a) of 
the band-heads. Once a band-head stands, an oscillation 
band will grow up from the band-head simply via the fluc- 
tuation of basic pairs. For examples, for L — 2 states, the 
dominant \a) of the band- heads of the four lowest oscil- 
lation bands are |(1) 2 P(°>), |(2) 1 p(°)), |(-2) 1 (1) 4 J P (0) ), 
and \(3) l (-l) l P^) with K a = 2,4,8, and 10, respec- 
tively. 

For L = 3 states, the dominant \a) of the band-heads of 
the three lowest bands are |(l) 3 p(°)), |(2) 1 (l) 1 p(°)), and 
|(3) 1 P^ ^), with K a — 3, 5, and 9, respectively. Since the 
p-, d-, and f-wave appear successively, these band- heads 
exhibit 1-cluster, 2-cluster, and 3-cluster structures, re- 
spectively, as shown in Fig. 4. 



-k, i.e., change the components |(fci) mi (k2) m2 P^) to 



|(-/ci) mi (-fc 2 ) m2 P (i) ), and so on. 

' Z,i 



Therefore \P 



(-£) 



{^ ( z\)\ and E(- L ' z <^ = E( L ' Z -' l \ 



Yrast line 



N=100, V =1, B =0.025 



FIG. 5: Energies of the yrast states with L — to 10. 

Let us study the yrast states ^j^, each is the lowest 
one for a given L. The energies of them are plotted in 
Fig. 5, their wave functions are found as 



(L) 
1,1 



(L,I,1) 



(l)L P (j) 



(9) 



where Aj i is very small. When L is small, the fluc- 
tuation of basic pairs is small, and the yrast states are 
dominated by the j = component |(1) L P(°'). When L is 
larger, the weight of the |(l) L p( ') component becomes 
smaller. E.g., when L = 0,2,4, and 10, the weights 
of |(1) L P(°)) are 0.94, 0.84, 0.75, and 0.54, respectively. 
Evidently, the energy going up linearly in the yrast line 
in Fig. 5 is mainly due to the linear increase of the number 
of p-wave particles. 



N=20, L=0 




L=3 



FIG. 4: p2 for the band-heads of L = 3 states, 9i = 0. 

Furthermore, a — L state can be derived from the 
corresponding L state simply by changing every k to 



FIG. 6: The spectrum of the L — states with N = 20, 
Vo = 1, and Orange = 0.025. Refer to Fig.l 

When N = 20, all the above qualitative features re- 
main unchanged. Examples are given in Fig. 6 and 7 to 
be compared with Fig.l and 4. Nonetheless, the decrease 
of N implies that the particles have a less chance to meet 
each other, thus the particle correlation is expected to be 
weaker. Quantitatively, it was found that (i) The spac- 
ing of adjacent oscillation levels becomes smaller, it is 
now ~2.2 to replace the previous 3.15 (ii) The fluctu- 
ation becomes weaker. E.g., the weig hts of |P( J ')) of 
the state are 0.01, 0.95, and 0.03 for j = 1,2, and 
3, respectively, while these weights would be 0.16, 0.46, 
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0.18 
0.16 
0.14- 
0.12- 



N=20, L=3 




60 120 180 240 300 360 

FIG. 7: The same as Fig.4 but with N = 20 

and 0.28 if N = 100. (iii) When N becomes small, the 
geometric features would become explicit. E.g., for the 
3-cluster structure, the difference between the maximum 
and minimum of p2 is ~0.007 in Fig.4, but <~ 0.033 in 
Fig.7. 



The decrease of V Q or 6 r 



was found to cause an 



effect similar to the decrease of N, the spectra would 
remain qualitatively unchanged. Quantitatively, when 
V a is changed from 1 to 0.1, the spacing inside a band 
is changed from ~ 3.15 to <~ 2.14, and the fluctuation 
becomes much weaker as expected. 

In what follows we study the vortex states. For an 
arbitrary L a < N/2, the spectra of the ^!^ f i ~ L °' > and 

^Zi^ states are found to be identical 15 , except the former 



shifts upward as a whole by N — 2L , namely, 
E (N-L„,z,i) = E (L ,z,i) + N _ 2Lq 



(10) 



Furthermore, their p2 are found to be identical. 

A A 

Let us define an operator X so that the state X\a) is 
related to \a) by changing every fcj in \a) to — fcj + 1, 
i.e., <p ki (6) to <p- ki +i{0) = e ie <j>_ ki (6). We further found 
from the numerical data that 



La) 



(11) 



holds exactly. In fact, X causes a reversion of rota- 
tion of each particle plus a collective excitation. It does 
not cause any change in particle correlation, therefore p2 
remains exactly unchanged. Thus the L large states, 
including the vortex states L — N, can be known from 
the L small states. 

The underlying physics of this finding is the separabil- 
ity of the Hamiltonian (it is emphasized that the sepa- 
rability is exact as can be proved by using mathematical 
induction). Let 9 co a — J^fOi/N, which describes a col- 
lective rotation. Then H = — ^ \- Hi nt , where Hint 

describes the relative (internal) motions and does not de- 
pend on 6 coU . Accordingly, E^ z ^ = L 2 /N + E^f A , 
the former is for collective and the latter is for relative 
(internal) motions. The eigen-states can be thereby sep- 



arated as 'F 



}_„iL8 c 
2¥ 6 



Vint 



The feature of the 



internal states i>^ l f'^ has been studied in [19]. Where 
it was found that, for an arbitrary L Q 



, (N+L ,Z,i) _ , (L ,Z,i) 
Vint ~ Vint 



(12) 



With these in mind, eq.(10) and (11) can be derived as 
follows. 

From the separability 



^(N-L) _ i i(N-L)0 coll , {N-L,Z,i) 



(13) 



When X acts on a wave function with L, from the defini- 

A 

tion of X, L should be changed to — L and an additional 
factor Y\ e 103 — e lN 9co " should be added, thus 

3 



X^ 



(L) 



i(N-L)6 c 



/2tt 



Vint 



(14) 



Due to (12), the right hand sides of (13) and (14) are 
equal, thereby (11) is proved. 



Furthermore, since ip int ' ' 

? (N-L,Z,i) _ jp(-L,Z,i) 
^int 



(ipini Z '^)*i tnc internal 
- E \nf' l) - Therefore, 
(N - Lf/N - L 2 /N = N — 2L. 



energy E mt 

E (N-L,Z,i) _ E (L,z. 

This recovers eq.(10), the energy difference arises purely 
from the difference in collective rotation. 

If the particles are tightly confined on the ring, rapidly 
rotating state with a large L = JN — L a would exist, 
where J is an integer. Their spectra would remain the 
same but shift upward by J(JN — 2L a ) from the spec- 
trum with L 



L , while ^ { zf Lo 



Xj^f°\ where 



Xj changes each 4> ki to 4>-ki+j- Thus the rapidly rotat- 
ing states have the same internal structure as the corre- 
sponding lower states but have a much stronger collective 
rotation. 

When TV increases greatly while V or 9 range decreases 
accordingly, the qualitative behaviors remain unchanged. 
E.g., when N = 10000 and V = 0.01 (0 range remains 
unchanged), the spectrum and the wave functions arc 
found to be nearly the same as the case N — 100 and 
V = 1, except that the spectrum has shifted upward 
nearly as a whole by 3939. This is again a signal that, 
for weak interaction and for the ground states, the mean- 
field theory is a good approximation. 

It is noted that the confinement by a ring is quite differ- 
ent from a 2-dimcnsional harmonic trap. In the latter, 
the energy of a particle in the lowest Landau levels is 
proportional to its angular momentum k. However, for 
the rings, it is proportional to k 2 . Consequently, higher 
partial waves are seriously suppressed and the p-wave ex- 
citation becomes dominant. For a harmonic trap it was 
found in [10,11] that d- and f-wave excitations are more 
important than the p-wave excitation when L is small. 
This situation does not appear in our case. 

When the zero-range interaction Vij — g5(0i — 9j) is 
adopted, The results are nearly the same with those from 
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the square-barrier interaction if the parameters are re- 
lated as g = 2V o range (in this choice both interactions 
have the same diagonal matrix elements). For an exam- 
ple, a comparison is made in Table III. The high simi- 
larity between the two sets of data imply that the above 
findings are also valid for zero-range interaction. 

TABLE III: Eigen-energies of the four lowest L — states 
for a system with N = 100 and with zero-range interaction 
Vij = 0.055(6*; — 8j) (the unit of energy is G as before). The 
weights of the j = components of these states are also 
listed. The corresponding results from square-barrier inter- 
action with V = 1, and Orange = 0.025 are given in the 
parentheses. 



(L,Z,i) 


E (L,Z,i) 


(cf ' z ' l) ) 2 


(0,1,1) 


39.0900 (39.0902) 


0.9370 (0.9371) 


(0,1,2) 


42.2982 (42.2983) 


0.0510 (0.0510) 


(0,7,3) 


45.4733 (45.4733) 


<0.02 (<0.02) 


(0,/7, 1) 


47.0076 (47.0074) 


0.8372 (0.8373) 



The numerical results from using zero-range interac- 
tion can be compared with the exact results from solving 
integral equations by Lieb and Liniger [16,17]. The vari- 
ables 7 and e(j) introduced in [16] are related to those of 
this paper as 7 = gir/N and e(7) = 4:ir 2 E/N 3 (the unit 
of E is G). However, this paper concerns mainly the case 
of weak interaction, say, g < 0.05, or 7 < 0.00157 (other- 
wise, the procedure of diagonalization would not be valid 
due to the cutoff of the space). Nonetheless, even 7 is 
as large as 0.5 (g — 15.9) the evolution of the ground 
state energy with N = 100 against 7 obtained via di- 
agonalization coincide, in the qualitative sense, with the 
exact results quite well . This is shown in Fig. 8 to be 
compared with Fig. 3 of [16], where 7 is ranged from 
to 10. In Fig. 8, the constraint 7 < e(j) is recovered. 
Furthermore, when 7 is small, e("f) against 7 appears as 
a straight line. 
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FIG. 8: e(7) = 4tt 2 E/N 3 against 7 = gir/N. N is given 
at 100 and E is the ground state energy calculated from the 
diagonalization in the unit G. 



In summary, a detailed analysis based on the numerical 
data of TV— boson systems on a ring with weak interac- 
tion has been made. The main result is the discovery of 
the basic pairs, which exist extensively in all the excited 
states and dominates the low-lying spectra. The fluc- 
tuation of basic pairs provides a common mechanism of 
oscillation, the low-lying states are thereby classified into 
oscillation bands. Each band is characterized by having 
its specific additional excitation of a few particles. Since 
the mechanism of oscillation is common, the level spac- 
ings of different bands are nearly equal in a spectrum. 

To divide the motion into being collective and relative 
provides a better understanding to the relation between 
the higher and lower states. The very high vortex states 
with L rj N can be understood from the correspond- 
ing low-lying states because they have exactly the same 
internal states. 

The particle correlation has been intuitively studied, 
particle densities are found to be in general non-uniform, 
bose-clustcring originating from the symmetrization of 
wave functions is found, which leads to the appearance 
of one, two, and three clusters. This phenomenon would 
become explicit and might be observed if N is small. 
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